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ABSTRACT 

We report a coordinated multi-band photometry of the RS Oph 2006 outburst and highlight the emission 
line free y-band photometry that shows a mid-plateau phase at y ~ 10.2 mag from day 40 to day 75 after the 
discovery followed by a sharp drop of the final decline. Such mid-plateau phases are observed in other two 
recurrent novae, U Sco and CI Aql, and are interpreted as a bright disk irradiated by the white dwarf. We have 
calculated theoretical light curves based on the optically thick wind theory and have reproduced the observed 
light curves including the mid-plateau phase and the final sharp decline. This final decline is identified with 
the end of steady hydrogen shell-burning, which turned out the day ~ 80. This turnoff date is consistent 
with the end of a supersoft X-ray phase observed with Swift. Our model suggests a white dwarf mass of 
1 .35 ± 0.01 M©, which indicates that RS Oph is a progenitor of Type la supernovae. We strongly recommend 
the y-filter observation of novae to detect both the presence of a disk and the hydrogen burning turn-off. 

Subject headings: binaries: close — binaries: symbiotic — novae, cataclysmic variables — stars: individual 
(RS Ophiuchi) — supernovae: general — white dwarfs 



1 . INTRODUCTION AND SUMMARY 

RS Oph is one of the well-observed recurrent novae and is 
sugges ted to be a progenitor s ystem of Type la supernovae 
(e.g., [Hachisu & Kato 2001b). It has under gone its sixth 
record ed outburst on 2006 February 12 UT (iNarumi et al.l 
2006). The five previous recorded outbursts occurred in 1898, 
1933, 1958, 1967, and 1985. These short 10-20 yr recur- 
rence periods indicate that the white dwarf (WD) mass is very 
close to the Chandrasekhar mass and that the mass accretion 
ra te is as large as M am ~ 1 x 10~ 7 Mq yr" 1 (see, e.g., Fig. 2 
of lHac hisu & Kato 2001b). If the WD mass increases after 
each outburst, RS Oph will eventually explode as a Type la 
supernova (e.g..lNomotolll982t lHachisu et al~ll 19961 1 1999a; 
lHachisu & Kato 1 1200 lbl) . Therefore, it is important to esti- 
mate the WD mass and the accreted mass left on the WD after 
the outburst. 

It is well known that the X-ray turnoff time is a good in - 
dicator of the WD mass (e.g., Hachisu & Kato 120051 12006a). 
When the hydrogen shell-burning a top the WD ex tinguishes, 
a supersoft X-ray phase ends (e.g. JKrautter et al7lll996l) . In 
a visual light curve, however, this turnoff is not clear because 
many strong emission lines contribute to it. To avoid such 
contamination to the continuum flux, we have observed RS 
Oph with the Stromgren y-band filter. The y-filter is an in- 
termediate bandpass filter designed to cut the strong emission 
lines in the wide V bandpass filter, so that its light curve rep- 
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TABLE 1 
List of Observers 



name of 
observer 


location 
in Japan 


telescope 
aperture 


observed 
bands 


No. of obs. 
nights (>■) 


Kiyota 


Tsukuba 


25cm 


B,V,y,R c ,I c 


24 


Kubotera 


Odawara 


16cm 


B,V,y,R c 


g 


Maehara 


Kawaguchi 


20/25cm 


B.V.y,R c ,I c 


19 


Nakajima 


Kumano 


25cm 


B,V,y,Rc,k 


55 


OKU a 


Kashiwara 


51cm 


V,y 


25 



a Osaka Kyoiku University team 



resents the continuum flux of novae (e.g., lHachisu & Kato I 
2006b). We have further modeled the light curve of RS Oph 
and have determined the WD mass by fitting our modeled 
light curve with the observation. 

Section |2] presents our multi-band photometry of the RS 
Oph 2006 outburst. The light curve fitting of the observation 
with our numerical model are presented in |5] Discussion 
follows in 50] 

2. OBSERVATION 

Optical observatio ns were started just after the discovery 
of the 2006 outburst (Narumi et al. 2006). Each observer and 
their observational details are listed in Tabled We have put a 
special emphasis on the Stromgren y-filter to avoid contami- 
nation by the strong emission lines. These y-filters were made 
by Custom Scientific Inc. 8 and distributed to each observer by 
one of the authors (M. Kato). Kiyota, Kubotera, Maehara, and 
Nakajima (VSOLJ members) started observation on February 
1 3 and obtained 65 nights data for y-magnitudes (from Febru- 
ary 17 to July 27). Osaka Kyoiku University (OKU) team ob- 
tained V and y magnitudes of 25 nights starting from February 
17 (until July 14). The magnitudes of this object were mea- 
sured by using the local standard star, TYC2 5094.92.1 (Kiy- 
ota) or TYC2 5094.283.1 (the other observers). We adapted 

8 http://www.customscientific.com/ 
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TABLE 2 
Parameters of RS Oph" 



parameter 


symbol 


25% 


50% 


100% 






efficiency 


efficiency 


efficiency 


WD mass 


M W D 


1.35M 


1.35M Q 


1.35M 


irradiation 


The 


0.25 


0.50 


1.0 


disk size 


a 


0.24(33 Rq) 


0.34 (47 R Q ) 


0.48 (66 Rq) 


size of RG 


7 


0.24 (25 Rq) 


0.34 (35 Rq) 


0.47 (48 Rq) 


distance 


d 


0.9 kpc 


1.3 kpc 


1.7 kpc 



a ^RG = O.7A/0, the inclination angle of the binary i = 33°, the separation a = 
316.5 R , R* = 138.3 « , fi 2 * = 102.5 R , /3 = 0.05, i/ = 2.0, T RG = 3550 K, the 
hydrogen burning turnoff date of day 83, and E(B — V) = 0.73 are common among all 
models. 



the brightness and color of (y = V = 9.57, B-V = 0.56) for 
TYC2 5094.92.1 and (y = V = 9.35, B-V = 1.23) for TYC2 
5094.283.1 from Tycho2 catalog. 

The y-magnitudes are plotted in Figure ^ together with I c - 
and V-magnitudes. We have also added visual magnitudes of 
the 1985 outburst from the American Association of Variable 
Star Observers (AAVSO) for comparison. Our y-magnitudes 
show very small scatter and follows the bottom of the 1985 
visual magnitude. The essential feature of the light curve is 
very similar to the previous outbursts. 

The y light curve, however, clearly shows a plateau phase 
from day 40 to day 75 and the sharp final decline starting from 
day 75. Such mid-plateau phases are also observed in two 
other recurre nt novae, U Sco (e.g.JHachisu et al. 2000) and 
CI Aal (e. g ..lHachisu&Katoir2001atlHachisuetal.ll2003t 
Hachisu & Kato 2003a). These authors interpreted the mid- 
plateau phase as a bright disk irradiated by the hydrogen- 
burning WD and a sharp start of the final decline as the epoch 
when the hydrogen shell-burning ends. 

3. LIGHT CURVE MODEL 

Our binary model, essentially the same as that in 
lHachisu & Kato I feOOlbl) except for the free-free emission 
model (see below) in the early phase of the outburst, consists 
of a red giant (RG) star, which is not filling its Roche lobe, a 
white dwarf (WD), and a disk around the WD. A circular orbit 
with the ephemeris given by Fekel et al. ( 2000) is assumed. 

3.1. Photospheric evolution of the white dwarf 

After a thermonuclear runaway sets in on a mass-accreting 
WD, its photosphere expands greatly and an optically thick 
wind massloss begins. The decay phase of novae can 
be followed by a sequ ence of steady state solutions (e.g., 
iKato & Hachisu if 1994) . After the optically thick winds stop, 
the envelope settles into a hydrostatic equilibrium where its 
mass is decreasing by nuclear burning. When the nuclear 
burning decays, the WD enters a cooling phase, in which 
the luminosity is supplied with heat flow from the ash of hy- 
drogen burning. We have followed nova evolutio n, using the 
same m ethod and numerical techniques as in lKato & Hachisu 
(1994). 

In our nova light curve model, we assume that free-free 
emission of the optically thin ejecta dominates the continuum 
flux in the early phase of RS Oph outbursts as in many classi- 
cal novae (e.g., Gallagher & Nev 1976). This is the main and 
most important difference from previous Hachisu & Kato's 
(2001b) model, in which the blackbody emission is assumed. 

The free-free emission of optically thin ejecta is estimated 
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FIG. 1. — Three (y, V, and I c ) bands light curves for the 2006 outburst 
of RS Oph. The AAVSO visual magnitudes of the previous 1985 outburst 
are added for comparison. We find a plateau phase from day 40 to day 75 
and a sharp final decline from day 75. Filled circles: /(-magnitudes. Filled 
triangles: y-magnitudes. Open diamond: V-magnitudes observed by OKU. 
Open squares: V-magnitudes observed by Kiyota, Kubotera, Maehara, and 
Nakajima (VSOLJ members). Small dots: visual magnitudes of the 1985 
outburst taken from AAVSO. 
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FIG. 2. — Calculated light curves for free-free emission during the op- 
tically thick wind phase together with observational y- (open circles) and 
/^-magnitudes (open squares) of the 2006 outbu rst. We also add observa- 
tional infrared K- (open triangles: taken from Evans et al, 1988) and vi- 
sual magnitudes (small dots: taken from AAVSO) of the 1985 outburst. For 
the y-magnitudes, we show six light curves for different WD masses, i.e., 
M W d = 1.3, 1.33, 1.34, 1.35, 1.36, and 1.37 M Q from right to left. For the 
I c - and K- magnitudes, we show two light curves for two WD masses, i.e., 
Mwd = 1.3 and 1.35 Mq from right to left. Here we assume X = 0.5 and 
Z = 0.02 for the chemical composition of the envelope. 
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FIG. 3. — Same as Fig. [2] but for a combination of our free-free plus disk 
irradiation models. The plateau phase is reproduced with a relatively large 
disk, a = 0.34, for the 50% iraddiation efficiency. The various parameters are 
summarized in Table|5] Thin solid lines denote the free-free emission bright- 
ness (labeled by "free-free") or the blackbody emission luminosity (labeled 
by "blackbody"). Thick solid lines are the total of the free-free plus black- 
body. Dash-dotted lines are the blackbody luminosity only from the WD and 
the RG with irradiation. Several V-magnitudes in the very early phase of 
the 2006 outburst (Open triangles) are added. We cannot reproduce the mid- 
plateau phase without a disk. The final sharp decline corresponds to the end 
of the hydrogen shell-burning. 



by equation (9) in Hachisu & Kato (2006b). The calculated 
free-free light curves are shown in Figure |2 The decline 
rate of the light curve, i.e., the evolutionary speed depends 
very sensitively on the WD m ass if its ma s s is very close to 
the Chandrasekhar mass (e.g., iKato II 19991 lHachisu & Katol 
2006a b). This is because the WD radius is very sensitive to 
the increase in mass near the Chandrasekhar mass. 

The timescale also depends weakly on the chemical com- 
position of envelopes. Hydrogen depletion is expected. This 
is because, just after the unstable nuclear burning sets in, con- 
vection widely develops and mixes processed helium with un- 
burned hydrogen. This mixing reduces the hydrogen content 
by 10%-20% for massive WDs like in RS Oph. The CNO 
abundance is not enriched in the recurrent novae, so we adopt 
the hydrogen content of X = 0.50 and the solar metallicity of 
Z = 0.02. 

We add ed infrared .^ -magnitudes of the 1985 outburst ob- 
served bv IE vans et al. I d!988h . Very little dependence of the 
light curve shape on the wavelength is a characteristic feature 
in the free-free emission light curves. However, the free-free 
light curve cannot reproduce the mid-plateau phase, so we in- 
troduce an irradiated disk in the next subsection. 

3.2. The irradiated disk and companion 
We assume an axi-symmetric disk with a size of 

Rdisk = aR*, (1) 

and a thickness of 

h = (3R dmk (®/R dmk y , (2) 
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FIG. 4. — The absorption-distance relation is plotted for the disk irradiation 
model and the companion model. Lines with open squares, open circles, and 
open triangles denote the disk irradiation model with the efficiency of 25%, 
50%, and 100%, respectively. Each symbol corresponds to the disk size, 
i.e., a = 0.64, 0.48, 0.32, and 0.24 from top to bottom. Lines with filled 
squares, filled circles, filled triangles, and filled stars denote the companion 
star model with the radius of 7 = 0.25, 0.30, 0.40, and 0.50, respectively. 
Each symbol corresponds to the photospheric temperature of the companion 
star, i.e., T Ra = 3200, 3300, 3400, 3500, 3600, 3700, 3800, and 3900 K from 
top to bottom. The absorption of Ay = 3.1E(B -V) = 2.26 is indicated by a 
dashed line 1 Sniiders 1987). Three models in Table l2lare indicated by large 
open stars. 



where /?disk is the outer edge of the disk, R* is the effective 
radius of the inner critical Roche lobe for the WD compo- 
nent, h is the height of the surface from the equatorial plane, 
tt! is the distance from the symmetry axis, and v is the power 
index which describes flaring-up of the disk e dge (see, e.g., 
iSchandl et al~IH997t lHachisu & Kato Il2001bl) . We also as- 
sume a companion red giant star with a radius of 



Arc = 7^2, 



(3) 



where is the effective radius of the inner critical Roche 
lobe for the red giant component, its mass of Mrg = 
0.7 Mm, and the inclination a n gle of the binary, i = 33° (e.g., 
IDobrzvcka & Kenvonlfl994l) . iDobrzvcka et al. I i ll 9961) sug- 
gested 7 ~ 0.4 for the distance of 1.5 kpc. 

The disk surface absorbs UV and supersoft X-ray photons 
from the WD and emits a part of it as a thermal emission 
with a lower temperature than that of the WD photosphere. 
The resultant light curve depends mainly on the disk size (a), 
the efficiency of irradiation (ry e ff = radiated energy/absorbed 
energy), and also the radius of the RG, but depends very 
weakly on the other two parameters of v and (3. Here, we 
assume v = 2 and (3 = 0.05. The dependence on these pa- 
ramet ers was widely discussed in the previous papers (e.g., 
lHachisu & Kato 1200 1 bl 12003 ^ . The irradiation efficiency of 
the RG hardly affects the light curve as shown in Figure|3] 

We have obtained three best fit models of the 2006 outburst 
in Table|3 The calculated light curves for the 50% efficiency 
are plotted in Figure [3] These models reproduce the mid- 
plateau phase and the sharp final decline identified as the end 
of hydrogen shell-burning. The turnoff date of day 83 in our 
1.35 M Q WD model is very consisting with t he supersoft X- 
ray tu rnoff on day ~ 90 observed with Swift dOsborne et al. I 

3.3. Distance 

The distance is obtained from the light curve fitting both 
at the plateau phase and at the post-outburst minimum phase 
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as shown in Figures [2] and [3] The disk luminosity depends 
mostly on the disk size, a, and the irradiation efficiency, ry e ff- 
We have changed these two parameters and calculated the 
brightness at the mid-plateau phase. Fitting the calculated 
brightness with the observation, we obtain the apparent dis- 
tance moduli, i.e., (m—M) y = (m—M)y and (m—M)j. Then we 
calculate the absorption from 

_ A V -A l _ (m-M)y-(m-M) l 

Av ~ 0.518 ~ 0.518 { ) 

where we use A/ = 0.482A y (e.g., Rie ke & Lebofskvlll985l) . 
Once Ay is obtained, the distance is calculated from log(ii) = 
((m — M)y— A v — 5)/5. Thus, we obtained the absorption- 
distance relation for the irradiated disk as shown in Figure 
We further r estrict the distance with the observed absorption 
of Av ~ 2.3 (Sniiders1 ll987l) . The same method is applied to 
the companion star, in which we have changed the companion 
size, 7, and the effective temperature, 7rg- The absorption- 
distance relation for the companion is also plotted. 

The largest ambiguity of our model is the irradiation effi- 
ciency of the disk. Here, the distances of 0.9, 1.3, and 1.7 kpc 
are derived for the three different assumed efficiency, i.e., 
25%, 50%, and 100%, respectively. The actual efficiency is 
somewhere between 50% and 100%, so we have a reasonable 
distance of 1 .3 - 1 .7 kpc. 

4. DISCUSSION 

The decline rate of free-free light curve and the hy- 
drogen burning turn-off date depend not only on the WD 
mass but also on the che mical composition of the envelope 
tHachisu & Kato l2006bt) . We have examined two other cases 
of hydrogen content, X = 0.70 and X = 0.35, and found 



that the best fit models are obtained for Mwd = 1-36 and 
1.34 M Q , respectively. So we conclude that the WD mass 
is 1.35 ±0.01 M Q . 

There are still debates on the distance to RS Oph. In the 
previous outburst, Hiellming et al. ( 1986) estimated the dis- 
tance to be 1 .6 kpc from H I absorption-line measurements. 
ISniiders I dl987l) also obtained the distance of 1 .6 kpc assum- 
ing the UV peak flux is equal to the Eddington luminosity. 
Harrison et al. ( 1993J) calculated a distance of 12 90 pc from 
the K-band luminosity. lHachisu & Kato I (|2p01b) obtained a 
smaller distance of 0.6 kpc from the comparison of observed 
and theoretical UV fluxes integrated for the wavelength re- 
gion of 91 1-3250 A . This shorter distance is caused by their 
blackbody assumption, because the flux is much larg er than 
the bla ckbody flux in this wavelength region. O'Brie n et al. I 
(2006) estimated the distance of 1.6 kpc from VLBA map- 
ping observation wi th an expansion vel ocity indicated from 
emission line width. Mo nnier et al.H2006l) estimated a shorter 
distance of < 540 pc assuming that the IR interferometry size 
corresponds to the binary separation. If we assume this cor- 
responds to a circumbinary disk, however, a much larger dis- 
tance is obtained. Therefore, our new value of 1 .3 - 1 .7 kpc is 
consistent with other estimates. 
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